Hybrid circuit cavity quantum electrodynamics with a micromechanical resonator 



J.-M. Pirkkalainen, S. U. ChoQ Jian Li, G. S. Paraoanu, P. J. Hakonen, and M. A. Sillanpaifj] 

O. V. Lounasmaa Laboratory, Low Temperature Laboratory, 
Aalto University, P.O. Box 15100, FI-00076 Aalto, Finland. 

Hybrid quantum systems with inherently distinct degrees of freedom play a key 
role in many physical phenomena. A strong coupling can make the constituents loose 
their individual character and form entangled states. The properties of these collective 
excitations, such as polaritons of light and phonons in semiconductors, can combine 
the benefits of each subsystem. In the emerging field of quantum information con- 
trol, a promising direction is provided by the combination between long-lived atomic 
states^! a nd the accessible electrical degrees of freedom in superconducting cavities 
and qubitPS. Here we demonstrate the possibility to integrate circuit cavity quan- 
tum electrodynamics^^ with phonons. Besides coupling to a microwave cavity, our 
superconducting transmon qubit interacts with a resonant phonon mode in a microme- 
chanical resonator, allowing the combination of long lifetime, strong tunable coupling, 
and ease of access. We measure the phonon Stark shift, as well as the splitting of the 
transmon qubit 7 spectral line into motional sidebands representing transitions between 
electromechanical polaritons formed by phonons and the qubit. In the time domain, 
we observe coherent sideband Rabi oscillations between the qubit states and phonons. 
This advance may allow for storage of quantum information in long-lived phonon states, 
and for investigations of strongly coupled quantum systems near the classical limit. 



Superconducting quantum bits based on Josephson 
junctions^ have offered an unparalleled testing ground for 
quantum mechanics in relatively large systems. At the 
same time, Josephson devices represent a promising im- 
plementation for quantum information processing. Basic 
quantum algorithms have indeed been recently demon- 
strated with both phase^ as well as transmoiPHED qubits. 
The latter operate in the circuit cavity quantum elec- 
trodynamics (QED) architechture, in which the qubits 
couple to an on-chip^ or 3-dimensional microwave cav- 
ity resonator^. The circuit QED setup, which enables 
coupling of qubits and non-destructive measurements of 
quantum states, can be regarded as the most feasible 
platform. 

The forthcoming challenges include the storage and 
retrieval of qubit states in a long-lived quantum mem- 
ory. Hybrid quantum systems, where macroscopic qubits 
are coupled to spin ensembles, are showing promise due 
to their long lifetime^EI], hut with the drawback of a 
difficult access and small coupling at the level of a sin- 
gle atomic degree of freedom. Another important chal- 
lenge is the quantum communication^^ between spa- 
tially separated superconducting qubits. 

Micromechanical resonators were brought to the quan- 
tum regime of their motion only very recently^ 17 ! 18 !, and 
their coupling to Josephson junction qubit S^HHJ has high 
potential to achieve the aforementioned goals. The first 
experiment of this kind demonstrated the interaction be- 
tween a charge qubit and a beam resonating at the fre- 
quency u) m /2ir ~ 60 MrkPH This experiment, using 
a read-out via the micromechanical device, did not al- 
low for time-resolved measurements, a pre-requisite for 
coherent state transfer. The resonant coupling was re- 
cently demonstrated with a microwave-regime piezoelec- 



tric device^, having w m /27r ~ 6.2 GHz, but the energy 
decay times were unfortunately very short due to high 
mechanical frequency and special materials. 

Our approach to improve on these issues is to utilize a 
membrane-type micromechanical resonator embedded in 
a complete circuit QED device. Our aim is also to further 
the understanding of the physics of slowly moving bodies 
close to the quantum limit, as well as the emergence of 
ultrastrong coupling of a micromechanical phonon degree 
of freedom, representing matter, to microwave light. 

On the conceptual level, our device is analogous to 
an optical cavity QED system where a two-level atom is 
coupled to two cavities having different frequencies (see 
Fig. la). We use a superconducting transmon qubit 7 , 
with the capacitance C 9 , single-electron charging energy 
E c ~ e 2 /2C 9 and Josephson energy Ej as an artifi- 
cial two-level system. It is coupled to a highly detuned 
phononic cavity formed by a suspended aluminum mem- 
brane (Fig. Id) with (jj m /(2w) ~ 72 MHz by means of 
a position-dependent capacitance C g {x), as well as to a 
nearly resonant photonic cavity realized as an on-chip 
coplanar waveguide microwave resonator with frequency 

U) c . 

The equivalent electrical circuit (Fig. lb) for the two- 
resonator circuit QED system allows writing the Hamil- 
tonian 

H =H q + Hu) c (a f a + 1/2) + hu m (Pb + 1/2) 

+ hg c (n - n)(a f + a) + hg m (n - h)(tf + 6), (1) 

where H q is the qubit Hamiltonian (see Supplementary) , 
a and b are the annihilation operators of the microwave 
cavity and of the mechanical resonator, respectively, and 
n is the number operator of Cooper pairs on the transmon 



2 



C Q 

► c coplanar waveguide cavity 

II— HE 




FIG. 1: Hybrid QED setup, a, Illustration of a cav- 
ity quantum electrodynamics (QED) system consisting of 
an atom (green) coupled to two cavities. Blue color cor- 
responds to the (photonic) microwave resonator and red to 
the (phononic) lower-frequency mechanical resonator, b, Cir- 
cuit schematics of an analogous electromechanical tripartite 
system, measured via a microwave cavity with a probe tone 
at frequency ui c . c, Optical image of the chip. Shown is 
the quarter-wavelength microwave cavity, the superconduct- 
ing transmon qubit playing the role of the atom, the mechan- 
ical resonator with its external control, a flux control line (or- 
ange), and the ground plane (magenta), d, Scanning electron 
micrograph showing the 5 fim long and 4 /im wide bridge-type 
mechanical resonator (dashed box) suspended above the qubit 
island. 



island. The qubit-mechanics interaction is given by the 
electromechanical coupling energy Hg m — x zp Vd c -^-^. 
Here, x zp is defined as the mechanical zero-point mo- 
tion, Vdc is a constant voltage applied to the mechanical 
resonator, and no — C g Vd c /2e is the dimensionless gate 



charge. The coupling between the qubit and the cavity, 
characterized by the energy hg c , allows for a dispersive 
measurement of the state of the qubit using the state- 
dependent pull on the cavity^. 

We use the three lowest energy eigenstates of the trans- 
mon qubit, labeled as the ground state \g), and the first 
and second excited states |e) and |/), respectively. The 
corresponding transition frequencies, tunable in the GHz 
regime by flux bias $dcj ar e ^> s - c and w e _f. Somewhat 
different from the standard concept of a qubit, which usu- 
ally involves the ground state, we treat |e) and |/) as the 
two-level system. This is beneficial, because the phonon 
coupling grows for higher levels of the qubit. 

The coupling of the qubit to the phonons in the me- 
chanical resonator in Eq. (1) is described by a generic 
light-matter interaction. The quantum states of this cou- 
pled system are dressed states, which are combinations 
of the qubit and of the mechanical resonator Fock states 
|iV m ) with phonon number N m . For high N m the ef- 
fective coupling g m \/N m + 1 becomes comparable to the 
qubit frequency, and the states attain a highly entangled 
character as electromechanical polaritons. Similarly, the 
qubit transition frequencies are expected to depend on 
the number of quanta (ac Stark shift). As a phonon- 
induced effect, it was measured indirectly through the 
mechanics^!. One can calculate the phonon Stark shift 
of the frequency w c _f of the |e)-|/) transition going be- 
yond the linear regime: 

Acj c _ f /27r = --y- cos(27m ) [J Q (2im x ) - 1] , (2) 

where the argument of the Bessel function is the mo- 



tional gate charge amplitude 



and 



the charge dispersion of the transmon qubit is e e -f — 

m _i_ 3 

Ec^p (fife) 2 3 exp (-^Ej/Ec) , with m = 
3 for this transition. The shift grows first linearly 
with the phonon number, according to Aw c _f ~ 

^tt 2 cos(2^n )(|^) 2 iV m . 

The device (see Fig. lc,d) was fabricated out of alu- 
minum on a sapphire substrate in a process comprising 
three layers of electron-beam lithography. In order to 
suspend the mechanical resonator, we used organic resist 
as a sacrificial layeJ^ (see Supplementary) . The quarter- 
wavelength microwave cavity has a resonant frequency of 
ui c /2tt = 4.84 GHz, coupling to the qubit capacitively 
at the rate g c /2Tr = 100 MHz. The cavity is coupled to 
the measurement circuitry by a large input/output ca- 
pacitor resulting in the external lincwidth ~/e/2t: ~ 15 
MHz. The device is mounted in a dilution refrigerator 
at a temperature of T = 25 mK. Although we did not 
access the thermal phonon number, we expect the mem- 
brane to thermalize^l down to the thermal occupancy 
= k B T/hv m < 10. 

We first identified the lowest mechanical eigenmode of 
the membrane. The electromechanical interaction was 
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FIG. 2: Phonon Stark shift, a, The driven mechanical resonance appears as enhanced absorption in a single-tone cavity 
spectroscopy, b, In two-tone spectroscopy (flux bias $dc/$o = 0.368, and no = 0.5) the dip representing the transition \e)- 
\f) (between the first and second excited states) of the transmon qubit, is observed to red-shift when the mechanical drive 
frequency is swept across resonance. Circles denote the phase of the cavity reflection, and the solid lines are Lorentzian fits. 
With resonant mechanical drive uj g /2n = u] rn /2n — 71.842 MHz, a mechanical motion is excited with an amplitude of 220 fm, 
about 80 times a; zp , and the phonon number N m ~6x 10 3 . In gate charge units, this equals a motional gate charge n x = 0.09. 
c, Shift of the qubit transition at flux bias &dc/&o = 0.370 as a function of driving frequency, with four drive amplitudes as 
labeled. Circles denote the center of the dip extracted from the data as in b, and the solid lines are the modeled response of a 
mechanical resonator with Q m — 5500. d, Highly non-linear regime of the phonon Stark shift with resonant drive at no = 0. 
Due to a single electron fluctuating on and off the island, also the transition line corresponding to no = 0.5 is visible as a mirror 
image. The black dashed line is a result of full numerical Floquet modeling, and the dashed white line is Eq. pj. 



switched on by applying a dc voltage V& c = 5 V, which 
permits an electromechanical coupling g m /2n = 2.5 
MHz. The mechanics was driven strongly by an addi- 
tional ac tone combined into a bias-T inside the cryostat 
as shown in Fig. lb. When the qubit is biased near the 
microwave cavity resonance, the cavity absorption in a 
single-tone measurement is altered when the mechanical 
mode is excited. This way, we obtain the mechanical 
peak at w m /27r = 71.842 MHz, as displayed in Figure 2a. 

By virtue of different frequency ranges of the qubit 
and the mechanical resonator, we cannot bring them 
on resonance by tuning the qubit frequency with a flux 
bias. A highly detuned regime, analogously to trapped 
ion d 1 * 26 * 2 -^, nonetheless allows for complete quantum con- 
trol. The pertinent phenomena are the Stark shift of the 
qubit transition frequency, and sideband transitions. 

A two-tone spectroscopy is used to investigate the 
qubit |e)-|/) transition for detecting the Stark shift. A 
standard method to probe this transition is to start by 
applying resonant microwaves in order to excite the qubit 
from the ground state \g) to |e). However, it turns out 
that in our system a small but nonzero population exists 
on |e) even in equilibrium, thus an initial excitation is 
not necessary. The actual excitation microwave tone (fre- 
quency w cx t) is applied to the transmon flux coil, and the 
increased population of the level |/) is distinguished as a 
phase shift of the probe tone. A shift of Aw c _f/27r = —18 



MHz in Fig. 2b from the bare value ui e - { /2n = 2.81 GHz 
can be attributed to the phonon- induced Stark shift. 

We can extract the center of the \e)-\f) transition dip 
and vary the motion amplitude, as shown in Fig. 2c. A 
slight difference of the bare value to that in Fig. 2b is 
due to a flux offset of ~ 2 m$o- Towards increasing vi- 
bration amplitude, the shift becomes more pronounced, 
and thereafter assumes a highly nonlinear and oscillatory 
behavior when the motional gate charge reaches values 
of the order n x > 0.5. This rather extreme regime can be 
explored by hxing a resonant mechanical drive, as imple- 
mented in the measurement shown in Fig. 2d. This data 
shows the Bessel-type oscillatory Stark shift, in accord 
with Eq. 

To model the shift, we first note that the total time- 
dependent gate charge is comprised of of two compo- 
nents: An g = n x + C g V r ac /2e, where the second term is 
a cross talk from the ac voltage drive on the mechanics. 
On resonance, the latter contributes 8 % of n x . Since 
the two terms have opposite phases when the drive is 
below/above u) m , they sum up to make the Stark shift 
asymmetric, as clearly seen in Fig. 2c. We used the mo- 
tional gate charge n x at the frequency w m as a common 
fitting parameter. The analytic result in Eq. |2]) already 
gives an excellent agreement with the measurement, as 
shown in Fig. 2d. Beyond analytics, we made a complete 
numerical calculation using a Floquet modeling of the 
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full transmon Hamiltonian (see Supplementary), match- 
ing well the observed spectroscopy in Fig. 2c, d. 

The excellent match to the theory clearly vindicates 
the phonon origin of the Stark shift. On top of that, the 
highly non-linear shift can be seen as a manifestation of 
ultrastrong coupling from light to mattei!^22l. Indeed, 
the frequency shift due to coupling goes far beyond the 
linear regime, exceeding the bare mechanical frequency. 

We can use sideband transitions to transfer q uanta 
between the qubit and the mechanical resonato j 26 * 27 l 
More specifically, the transitions occur between the en- 
tangled dressed states. However, an intuitive appealing 
picture is to describe them by using uncoupled states 
\e,N m ), \ f 1 N' m ), see Fig. 3a. With the parameter choice 
in Fig. 3b, the transmon charge dispersion and effective 
coupling are large enough so that the basis states are en- 
tangled electromechanical polariton states. In Fig. 3c, 
the effective coupling is smaller such that the description 
in terms of separable qubit and phonon states is valid. 

Since the data of Fig. 3b were taken at the gate charge 
extreme offsets (no = 0, 0.5), we obtain only two-phonon 
processes due to symmetry reasons. However, offset drift 
during data taking rendered the first order sidebands vis- 
ible in Fig. 3b. As the Rabi frequency is relatively low 
in Fig. 3b, we can use an analytical result for the qubit 
population, obtaining an excellent agreement with the 
measured sidebands. Under proper phonon occupancy, 
we also can eliminate dephasing due to gate charge fluc- 
tuations, as well as the fluctuating quasiparticle. 

For investigating the energy exchange in the time do- 
main, we performed time-resolved measurements in the 
setup of Fig. 3c. Here, instead of a continuous microwave 
irradiation as in the two-tone spectroscopy of Figs. 2 and 
3, we pulsate the excitation microwave tone at varying 
widths, while monitoring the level |/) population. Dur- 
ing the pulse, the system is expected to coherently evolve 
at the Rabi frequency (for sideband k) between a pair 
of eigenstates as illustrated by arrows in Fig. 3a. Fig- 
ures 4a,b show the obtained Rabi oscillations on the bare 
qubit transition (k = 0) and on the first blue sideband 
(k = 1). The latter indicate that the qubit can swap 
the state while simultaneously adding or removing one 
quantum in the mechanical resonator, thereafter, con- 
verting the information into flying microwave photons in 
the non-destructive measurement process. We can map 
this energy exchange by varying the detuning of the ex- 
citation microwave (Fig. 4c), as well as by converting the 
data to frequency domain (Fig. 4d). The data shows, as 
expected, an increase of the Rabi frequency on detuning 
from a particular sideband. However, the same theory, 
which supposes non-overlapping of the peaks, cannot be 
used to relate the peak widths in Fig. 3c to the Rabi 
frequencies. 

The phonon-photon conversion demonstrated here 
forms a building block for the conversion of quantum in- 
formation between microwave light and mechanical mo- 
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FIG. 3: Motional sideband transitions, a, In a Stokes 
(red sideband) scattering event, the qubit exchanges (ab- 
sorbs/emits) a quantum of frequency w e — f + kuj m with the mi- 
crowave field, whereas the mechanics emits/absorbs k quanta. 
In an anti-Stokes (blue sideband) scattering event, the pro- 
cess is reversed, b, Cross-cut along the dashed vertical line 
from Fig. 2d at n x ~ 0.6. The solid black line is a fit from the 
analytical model presented in the Supplementary, with decay 
rate ^y/2n = 3 MHz, dephasing rate 7 v /27r = 6 MHz, and 
main peak Rabi frequency Ao/27r = 7 MHz. c, Spectroscopy 
with the motional gate charge corresponding to the dynamical 
sweet spot (n x = 0.4, see Supplementary), where the qubit is 
insensitive to offset charge fluctuations, $dc/$o = 0.326. The 
qubit was driven hard at Ao/27r = 31 MHz, making the main 
peak broad. 



tion. It also opens up the possibility of quantum mem- 
ory operation. While here the storage time is lim- 
ited by the thermal coherence time of the mechanics 
tt = Qm^m/N^ ~ 1 jits, we note the encouraging recent 
findings^ with Q m ~ 10 10 and tt of the order seconds. 
This kind of micromechanical device could further allow 
for quantum communication by converting phonons into 
flying photons using radiation pressure coupling, natu- 
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FIG. 4: Electromechanical Rabi oscillations, a, Oscillations of the phase shift of the probe tone, proportional to the 
population of the qubit third level |/}, measured near the bare qubit frequency k — 0. Traces are offset for clarity, b, Time- 
resolved evolution between phonons and qubit states, near the first blue sideband k = 1, with the resonant Rabi frequency 
Ai/2n ~ 9.4 MHz. The quantum states on the course of the evolution are marked, c, Qubit population oscillations mapped as 
a function of the qubit excitation frequency and Rabi pulse width. The traces in a and b are along the arrows. The decay time 
constant, related to total decoherence rate, is ~ 70 ns. d, Fourier transform of c showing the Rabi frequencies of each sideband 
growing as a function of excitation detuning. The lines are fits to the expression Afc(wext) = vA~I + ( u e-f — w ox t — kuj m ) 2 - 



rally suited for a membrane mirror. An increase of the 
electromechanical coupling by means of a narrower vac- 
uum gap and a smaller qubit capacitance could progress 
the interaction between light and true matter all the 
way into the single-phonon ultrastrong regime. The mo- 
tional sideband operations also facilitate the use of the 
full quantum toolbox developed in the field of trapped 

iQn jl | 26 | 27 | 

for engineering non-classical states of mechan- 
ical resonators, or long-distance entanglement of several 
slowly moving purely mechanical objects. 
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DEVICE FABRICATION 



The fabrication (Fig. SI I includes three layers of electron-beam lithography on a sapphire substrate. The first 
lithography patterns everything else except the mechanical resonator. Aluminum is deposited by shadow evaporation 
at thicknesses of 20 nm and 40 nm, with an oxidation in between in order to create the Josephson tunnel junctions. 

The sacrificial layer separating the bridge from the transmon island is defined with PMMA used as a negative 
resist. We start by spin-coating 3 % 950k PMMA, creating about 200 nm thick layer. Under high electron dose, here 
9 mC/cm 2 at 200 pA, PMMA cross-links, shrinks in thickness down to about 40 % and becomes insoluble. After the 
exposure, the non-exposed PMMA is washed away with acetone. 

For the third lithography, a 800 nm thick resist serves as the deposition mask for 300 nm Al defining the mechanical 
resonator. After lift-off, 90 min of isotropic O2 ashing at 1 Torr removes the PMMA and suspends the bridge. 
Although the ashing takes time, no appreciable change to qubit properties can be attributed to it. 

A highly tilted SEM micrograph (Fig. Sli>) reveals an undulating vacuum gap of about 40. ..100 nm between the 
bridge and the qubit island. Apparently, the PMMA thickness became non-uniform during the cross-linking. This 
makes it difficult to calculate d x C g from the geometry, however, the fitted value d x C g = 26 nF/m is in the proper 
range. 



Al 60 nm 

Sapphire 

negative PMMA 40...80 nm 





FIG. SI: Device fabrication, a, Steps (from top to bottom) of the process, b, Gracing angle scanning electron micrograph 
of the mechanical resonator (red) above the superconducting qubit island (green). 



BASIC CHARACTERIZATION 



With the electromechanical coupling set to zero by Vdc = 0, we first characterized the basic operation of the 
transmon-cavity circuit QED system. In Fig. S2 1 we show a single-tone spectroscopy, with the avoided crossings at 
roughly ^dc/^o — ±0.1, where the qubit frequency crosses the cavity. In fact, the crossing just barely takes place, 
since the maximum qubit frequency is 5.0 GHz at <&dc = 0. The wings, absent in the ideal case, around the crossing 
points are attributed to higher transitions in the coupled system. We expect them to be due to microwave leakage 
via the high gate capacitance. 
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Two-tone spectroscopy is the main tool used in the measurements described in the main text. It is the best way to 



reveal the un-dressed qubit spectral lines when far-detuned from the cavity. Fig. S2 d displays such characterization 
over a large span of flux biases and qubit frequencies. By overlaying the results from numerical modeling, the spectral 
lines can be distinguished as |g)-|e) and \e)-\f) transitions. We adjusted the gate voltage such that no = 0, however, 
especially the |e)-|/) line can be identified having split into two. The other branch (no = 1/2) corresponds to an extra 
quasiparticle on the island, excited likely by microwave leakage. 



|Sj(dB) 



arg(S 11 )(deg) 




0.32 0.34 0.36 0.38 
<& e /<& 



FIG. S2: Qubit spectroscopy, a, One-tone spectroscopy of the cavity resonance absorption, while the effective Josephson 
energy is tuned by flux bias, b, Two-tone spectroscopy, where the probe tone uu c is kept fixed near the bare cavity frequency 
(4.84 GHz), and qubit excitation frequency is swept. The dashed lines are the calculated transitions as labeled. 



The data of Fig. [S2] also allows for accurate fitting of the qubit parameters; we obtain Ej\ = 4.63 GHz, Ej2 = 6.43 
GHz, C q = 61 fF, Ec = 318 MHz, and the coupling energy to the microwave cavity g c ~ 100 MHz. 

The dc voltages applied to the mechanical resonator (acting also as a gate electrode to the qubit) make the 
background charge of the qubit more prone to jumps than without the dc voltage. We applied a slow digital feedback 
about once a minute to the gate in order to correct for the jumps. For this, we repeated a sequence consisting of 
measurement of the gate periodicity, followed by a small voltage offset, until the original offset was reached. 

The mechanical Q- value of the present sample Q m ~ 5500 is clearly lower than typically recorded (> 10 5 ) with 
aluminum beams^or membranes^ at cryogenic temperatures. We believe this is due to either clamping losses, or the 
oxygen ashing changing material properties. However, the NIST desigrPfor aluminum membrane, with Q m ~ 3 x 10 5 
could equally well integrate in the present setup, allowing for longer storage times of qubit states in the mechanics. 



DATA ANALYSIS 



The gate capacitance to the qubit, C g — 1.76 fF, is extracted from the gate charge no = C s Vd c /(2e) periodicity 
of the qubit transitions. Using the modeling of Stark shift asymmetry in Fig. 2c in main text, we can deduce the 
gate charge amplitude due to the direct driving. The driving voltage amplitude can then be calculated with the 
capacitance to the qubit. Therefore, the fitted gate charge amplitude due to the mechanical capacitance modulation 
translates into the derivative of the capacitance with respect to the distance, d x C g = 26 nF/m. This in turn gives 
the coupling between the qubit and the mechanical resonator to be g m /2Tr — 5.0 MHz with Vdc = 10 V which was 
the maximum voltage applied to the device. The obtained values of the driving voltage agree well with the estimated 
attenuation of the drive cable, although the latter has inaccuracies of several dB. 



Next, we discuss additional data on the phonon Stark shift. In figure S3 the gate charge is tuned to no = 0.25, at 
which point the qubit transition frequency is insensitive to quasiparticle jumps and only one transition line is observed. 
Equation (2) in main text, indeed, gives no dependence of the transition frequency on the modulation amplitude. 
Accurate numerics, plotted as a dashed line, however, displays some dependence owing to deviations from simple 
sinusoidal charge dispersion. Due to strong 1/f charge noise when off the gate charge sweet spot (no = 0;0.5), the 
transition peak as well as the phonon sidebands are nearly absent except at the dynamical sweet spots (see Fig. |S5[j) . 
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FIG. S3: Phonon Stark shift. Similar to Fig. 2d in main text, but the gate offset is set to the maximum slope, no = 0.25, of 
the modulation curve. The measured population of |/} is proportional to the color code. At this gate offset, the Stark shift is 
minimal as expected from both analytics and numerics. The dashed line is from full numerical Floquet modeling. 



THEORETICAL DETAILS 



Hamiltonian 



We are considering the tripartite system consisting of the transition qubit, microwave cavity, and the mechanical 
resonator, as displayed in Fig. |S4| 

The qubit is connected to the mechanical resonator via the gate capacitance C g {x) dependent on the displacement 
x of the mechanical resonator, which in the following can be either a classical variable or an operator. In the latter 
case, we have 



x zp (tf+b). (SI) 



The zero-point oscillation rms amplitude is 



2mw m 

The displacement gives rise to a motional gate charge defined as 

rix= dC M^L x . (S2) 
ax 2e 

The transmon qubit^ is schematically the same as the split Cooper-pair box, consisting of a superconducting loop 
interrupted by two Josephson junctions with Josephson energies Ej\ and £j2- For future use, we define the total 
Josephson energy Ej, and the difference Ej^. They are tunable by the flux <I> e through the superconducting loop: 

Ej(<f> e )=(Ej 1 +E J2 ) cos (tt^ 

} ?( P» 

Ej_ (* e ) = [En - E J2 ) sin tt^ 
For the moment, we linearize the dependence of the capacitance on displacement: 

CM = C g + ^lx = C g (l+ gx) , (S4) 



where we denoted the un-displaced gate capacitance by C g (x = 0) = C g . The capacitance in parallel with the 
junctions is denoted by C q , the sum of the capacitances of the individual junctions and of an interdigital capacitor. 



10 




FIG. S4: Circuit schematics of a transmon qubit (shaded box) coupled capacitively to a micromechanical resonator via a gate 
capacitance C g (x), as well as to a microwave cavity represented as an LC resonator (dashed box). 



The total capacitance of the qubit is the sum of C q , C g [x), and of C t which is the coupling capacitor to cavity. The 
values are supposed to be related as C g ,C t <C C q . The (single-electron) charging energy is 

e 2 e 2 
E ° = 2[C q + C g (x) + C t ] ~ 2C~ q " 

Let us derive the Hamiltonian for the entire tripartite system. The canonical coordinates tp and denote the phases 
of the qubit island, and of the cavity, respectively. The qubit island charge Q is related to the number of Cooper pairs 
n on the island as 

Q = 2en = C g (tp - V dc ) + C t (p -6) + C q tp . 

One writes down the Lagrangian which is a function of the generalized coordinates and their time derivatives. We 
use the subscript q for the qubit, c for the cavity, and m for the mechanical resonator. The Hamiltonian will consist 
of the Hamiltonians of the qubit H qi the cavity H c , the mechanical resonator H m , and of their couplings. We expand 

in the small parameter ^r-x, and obtain 

H =H q + H c + H m + H qc + H qm + H cm + H qcm + H c2m ■ (S5) 

The terms are 

H q = AE C [n - n ) 2 - Ej cos (tp) - Ej^ sin (ip) 
H c = hLu c (a) + a) 
H m =tm m {tf +b) 
H qc = hg c (n - n)(a) + a) 
H qm = hg m (n - n)(tf + b) 
H cm = frgcm(a^ + a)(6 f + b) 
H qcm = hg qcm n(a^ + a)(fe f + b) 
H c2m = %c2m(a t + a) 2 {tf + b) , 



(S6) 
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and the coupling energies read 



hg c = 
hg m = 



C c e /2hw 
~C, 



c 

, 2e 

X Zp ^dc Cg s Y 



x zp C'C c V dc j2huj, 



2C n 



c 



(S7) 



C2 



X Zp C g hUJ C 

%c2m — 



2C 



The direct cavity-mechanics linear interaction if cm is usually negligible, owing to different energy scales, 
The last term H C 2 m m Eqs. (S6) and (S7) is a radiation-pressure interaction from the mechanics to cavity. Although 



generally small, it can provide a medium for finding and characterization of the mechanical resonance via the "circuit 
optomechanics" methods^ by pumping the cavity up to photon numbers 3> 1. The second-last term H qcm also becomes 
in the dispersive limit approximately a radiation-pressure term™ however, it becomes irrelevant since a high photon 
number destroys the qubit operation. 

The relevant physics discussed here is thus due to the Hamiltonian of Eq. (1) in the main text, 



H — H a + H c + H„ 



(S8) 



which describes the qubit coupled to two harmonic oscillators. 

Let us repeat that the above results were written under the assumption of linearization of the gate capacitance 
with respect to displacement of the mechanical resonator, Eq. |S4|. This holds for small displacement such that the 
motional gate charge n x is much less that the period of 1. However, in the experiment we reach the regime where the 
linearization fails. The interesting physics in this regime is detailed starting from section [| 



Quantum treatment of the electromechanical coupling 

The qubit-mechanics interaction part in Eq. (1) in the main text is (setting h = 1) 

H = H q + Lu m (tfb + 1/2) + g m (n - n)(6 f + b). (S9) 

One can diagonalize the qubit, and rotate the coupling term into the qubit eigenbasis. Generally, the coupling will 
then have components in all three pseudospin directions. However, for this example we write it as 

H ~ -^Tz + u m {tfb + 1/2) + <? m (&t + b)a x , (S10) 

which is the Jaynes-Cummings Hamiltonian. The eigenstates are dressed states, which are entangled states of the 
qubit states, and of the harmonic oscillator Fock states \N m ) with N m phonons: 

|+, N m ) = cos (®^\ I/, N m ) + sin \e, N m + 1) 

\-,N m ) = - sin (^j |/, N m ) + cos (^j \e, N m + 1) . 

Here, 



The eigenenergies are 



, Q 2g my /N m + 1 

tan Bat = — , A = cj c _f — oj m \pl2) 



(N rn + 1) ± l -^AgUN m + 1) + A2 . (S13) 
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The degree of mixing is given by the angle &n which basically depends on the ratio of coupling to the qubit 
frequency, and on the number of phonons. By increasing the phonon number, therefore, the eigenstates can attain 
a highly mixed character. Another way to indicate the increased coupling with higher Fock states is that from 



Eq. (S12), ( S13 1 one sees that an effective coupling is 



.9 =9, 



which in the present case, given a maximum about N m ~2x 10 6 , would exceed the qubit frequency. 



The discussion starting from Eq. (S9| assumes the validity of the linearized interaction Eq. (S4|, however, this 
assumption does not hold any more when the motional gate charge becomes comparable to about n x > 0.1. A rough 
estimate for the mixing angle is then given by using for the effective coupling the amplitude for the transmon charge 
dispersion, see below Eq. (S25), viz., g* ~ e„ ~ 200 MHz. We then obtain that the dressing of the eigenstates remains 
quite moderate, at the level of ~ 10 % in the experiment. In the first approximation, the eigenstates are hence states 
with the qubit in either \e) or |/), and the mechanical resonator with a certain N m number of quanta. 

In order to treat effects beyond the linear regime, we will consider the coupling of the qubit to the mechanical 
resonator semiclassically, i.e., treating the latter as a classical field coupling to the qubit: 



n g(t) = n o + n x cos(ujgt) . 
This is the starting point for the rest of this Supplementary. 



(S14) 



Charge qubit approximation 



In the basis of number states, the general Hamiltonian of the split Cooper-pair box, or transmon, reads 

1 



Y,{ AE c (« - n g f \n) - -Ej (|n + 1) + \n - 1)) - % -Ej- (\n + 1) - \n- 1)) }(n| . 



(S15) 



Let us first consider the case of a charge qubit, which has Ec ^ Ej. Although the transmon qubit operates in the 
extreme limit Ec <C Ej, the present discussion provides, nonetheless, useful insight. For the moment, we ignore the 
asymmetry of the Josephson energies. 



Restricting to a subspace involving two states, Eq. (S15) can be written using the Pauli matrices. With the two 
lowest-energy states n = 0, 1, as well as n g between and —1/2, we obtain 



-2E c (l 



2n g )a z 



1 



-E,a 



(S16) 



The Hamiltonian Hqi has the eigenstates \g) and |e), the ground state and the first excited state, respectively. Of 
most interest for the present work is the next higher transition: 



H 



12 



-&E c n g a z 



-E in 

2 



J&X 



(S17) 



which has the eigenstates |e) and |/). Hence, the effective charging energy, and thus the charge dispersion of the 
energy, is enhanced by a factor of two by considering the first excited and second excited charge states. Generally, it 
grows by a factor of two for each higher transition in the charge qubit model. Notice that the gate charge periodicity 



in Eq. (S17) is offset by half a period with respect to the lowest transition, Eq. ( S16 1 



We swap the z and x spin indices in Eq. (S16l, for the purpose of following the standard that the spin is pointing 
along the z axis. The coupling of the driving field via charging energy is then perpendicular to the pseudospin 
direction. The two-charge state Hamiltonian describing the ground and first excited states g and e is: 



#01 



1 



Eja z -2E c (l + 2n g )a x 



(S18) 



with n g between and 



-1/2, see Fig. S5 



Beyond these values setting the limits of the linear regime (in the charge 
qubit case), the charging energy is periodic. We can approximate this dependence as sinusoidal, such that Eq. (S18) 
becomes 



Hi 



1 



01 



Ejcr z - E c [1 + cos(27TO g )] a x . 



(S19) 
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FIG. S5: Charging energy of the qubit. a, Charg ing energy for the two lowest charge states in Eq. ( S16 1 (black solid line), 
and a sinusoidal approximation of it (red, Eq. (S19l). The dashed lines depict the presence of an extra quasiparticle on the 
island. The y-axis is scaled with Ec and x-axis is the gate charge, b, Illustration of a "dynamical sweet spot" , where the qubit 
level spacing at a certain modulation amplitude is insensitive to changes in the offset charge no. The lines depicting the level 
spacings cross, independently of no, for the first time at about n x ~ 0.4. 



Using Eq. (S14), we obtain all harmonics by the Bessel expansion of the sinusoid of a sinusoid, while the dominant 



term is time-independent: 



-Era 



E, 



jl + cos(2Tm )J (2Trn x )^<7 x 



(S20) 



The difference between level spacings of Eq. (S20) and that of the uncoupled qubit, Eq. ( S19 1 with n g {t) = no, yields 
the Stark shift for the transition between the ground state g and the first excited state e: 



2E% 



Jg-e 



E 



J 



2cos(27rn )(Jo(27rn I ) - 1) + cos 2 (27m )( Jq (2im x ) - 1) 



(S21) 



Similar to Eq. (S19), we can approximate the higher transition Eq. (S17) 

Ej 



Hi 



-a z — 2E C [1 — cos(27rn ff )] a x 



(S22) 



and we obtain a Stark shift for the next higher pair of states e and /, enhanced by a factor of 4 over Eq. (S21 ): 

2cos(27rn )(Jo(27m ;r ) - 1) + cos 2 (27m )( Jq {2nn x ) - 1) 



Au e _ f 



8E^ 
Ej 



(S23) 



The results in Eqs. ( S21 ), ( S23 ) arise from a simplified model of treating the transmon as a charge qubit, although the 
parameter regime Ej/Eq of these two are in fact the opposite limits. However, the results are qualitatively correct, 
as seen in Fig. S6 This is mostly due to the fact that the charge dispersion for the higher transition (see section [| is 
relatively correctly predicted. 

In the experiment, we also observe the effect of a quasiparticle which intermittently tunnels in and out of the 
island. The presence of the quasiparticle shifts the charging energy by half a period, see Fig. |S5^ i. Therefore, in a 
time- averaged measurement, two gate charge values differing by one electron appear as superimposed on top of each 
other. This behavior is rather typically observed in circuits where the single-electron charging effects are noticeable^. 
As it is often the case in superconducting devices, the origin of the quasiparticles is somewhat unclear, but can be 
associated with radiation from higher-temperature parts of the experimental setup. 



Using the charge dispersion of the transmon qubit 

The charge dispersion of the transmon is derived in the limit Ej ^> E(p*. Here, for convenience we use numeric 
labels n = 0, 1, 2, ... for the eigenstates (not to be confused with number states used in sections [J |) . The eigenenergies 
are 

= w„(n ff = 1/4) - y cos(27m g ) . (S24) 
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The peak-to-peak amplitude is 



2 4 "+ 5 [2f Ej 



c 



exp 



(-y/8Ej/E C 



(S25) 



Since the higher levels have progressively stronger charge dispersion, a good approximation for the transition energies 
between two levels is given by tJ n -i <n — w n . 

Let us take the gate dependence as a slow modulation of the level energy in the eigenbasis, and, for simplicity, label 
cjq = u n (n g = 1/4). We use the previous notation for the Pauli matrices which operate now in the eigenbasis, 



H, 



q - -- [ujq + e 2 cos (2im g (t))} a z 



(S26) 



where n g (t) is given by Eq. (S14|. In the same way as starting from Eq. (S19) in the charge-qubit model in section 
[J we use a Bessel expansion to break this into multiple frequencies. However, in contrast to the simplistic model, we 
now include all frequencies: 



H q = -~S a z - ~ y^ff cos (fiuigt) <t z 



(S27) 



0=1 



where 



S =^o + — cos(27rn ) Jo(27rn x ) 

5/3 =(— l)' 9 ^ 2 e2 cos(27rr7, ) Jp(2irn x ), even /3 
5p =(-l) (/3+1)/2 e 2 sin(2™ )J^(2^), odd (3 



(S28) 



Spectroscopy by the excitation microwave 

First sidebands 



Let us consider as an introductory example a two-level Hamiltonian in its eigenbasis, coupled to two linearly 
polarized fields: 



H q = [S a + S cos(u g t)} a z - g cxt cos(uj cxt t)cr x . 



(S29) 



The (low-frequency) field with frequency uj g and Rabi frequency 6 will describe the (strong) drive due to the phonon 
field, and the perpendicular field having the frequency w cxt and Rabi frequency g cxt models the excitation tone in 
two-tone spectroscopy. The probe tone used to measure the qubit state via dispersive cavity readout is not considered; 
it simply gives a phase shift proportional to the population of the upper level. 

We move to a non-uniformly rotating frame with respect to the mechanical modulation by the transformation 



u = exp 



2oj, 



sin(w s i)<7 2 



We obtain 



H q = - —cr z - 22 cos [(w oxt + 2kuj g )t] a x + ]T A 2fe _i sin [w cxt + (2k - l)u> g t] a v . 



k= — oo 



k= — oo 



The Rabi frequencies are 



(S30) 



(S31) 



(S32) 
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Strong drive 



In reality the situation is more complex, since our two-level system (TLS) is coupled to a sum of fields having 
different frequencies, which are due to expanding the sine of sine (Eq. (|S27|), and also to the excitation tone: 



H « = -2 



S o + ^ S P cos(/3ujgt) 



8=1 



Oz - 5cxt cos(cJe X tJ)°V 



(S33) 



Here, one has to repeat the transformation u of Eq. ( S30 ) separately for each frequency component, yielding the 
overall transformation: 



U = exp 



2 



Equation (S33) becomes 



H q = --S a z - cos(w cxt t)(P + P*)o- x - cos(w oxt £)(P - P*)a y , 



(S34) 



(S35) 



where we used the shorthand notation 



P = 



exp 



3=1 1 9 



Let us also define = . With further algebra we get 



(S36) 



p* 



f 6 \ 

exp(— ipujgt) x J q ( — J exp(— 2iquj g t) x 



p— — 00 ^ ' q— — 00 ^ ' 

= |... + J-i(S'i) exp(ia; g t) + Jo{6[) + Ji(5[)cxp(-iuj g t) + J 2 (S[) cxp(-2iuj g t) + ... jx 
x j... + J-i{5' 2 ) exp(2iuj g t) + J (5' 2 ) + Ji(S' 2 ) exp(-2iw g i) + J 2 (5' 2 ) exp(-4iu; g t) + ...j) 



(S37) 



We neglect sidebands of sidebands, since they have a lower amplitude. Then the explicitly written two lines 
Eq. (S37) yield multiples of the mechanical frequency up to third order: 

P* ~ J (6[)J (S 2 ) - 2iJ 1 (6' 1 )J (5 2 )sm(u g t) - 2iJ 1 (6' 2 ) J (S[) sin(2u; 9 i) - 2iJ 1 (S' 3 )J (5' 2 )sm{3io g t) ' 

This result is applied in Eq. ( |S35 ), obtainin 



111 



(S38) 



H q = -~5 a z - A cos(w cxt t)a x - A x [sin(w oxt + u g t) - sin( 
+A 2 [sin(cj ext + 2uj g t) - sin(o; ext - 2u g t)] a y + ... 
The Rabi frequencies are now 



W eX t - Wgt)] <7y + 



(S39) 



A =. 9cxt Jo Jo (^) 

Al =9cxtJl (|) J ° (^) 

A2 (^) Jo (|) 

A 3 =5cxtJi(^-) Jo(^ 



(S40) 



A3 ^ (t g ) J ° (£) ■ 

The higher Rabi frequencies (k > 1) are suppressed from those in Eq. (S32), nearly comparable to that of the main 
peak (k = 0), due to the nested Bcssel function dependence. 
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Level populations 

One can estimate the populations of the qubit levels in the spectroscopy experiment discussed in section]] Let us 
consider a generic TLS described by a spin Hamiltonian, with a notation reminiscent of this discussion: 



H q = -~(T z -^-a x -—a v . (S41) 



This Hamiltonian is time-independent, however, in a rotating frame, a driven system will adopt similar form. For 
notational simplicity, let and 1 denote the ground and first excited states of the TLS. 
The Liouvillean master equation for the density matrix is: 

P = 4 [H, p] + C\p] , (S42) 
in 

where the Liouvillean with decay and pure dephasing is: 

£\P\ =~\ {?+°~P + P° + °~ - 2a-pa+) + ^ {a z pa z - p) . (S43) 
The steady-state population of the excited state in the Rotating- Wave Approximation becomes 

(A 2 X + A 2 y ) 

/°n = 1 - Poo = "To ^r-l , . , t^n s- , (S44) 

S 2 + 7^ ( A 2 + A 2) + 7 2 Qt 



Here, the total decoherence rate is 



7tot = ^ + 70 • (S45) 



Equation ( S44 ) is Lorentzian having the full width at half maximum 



^ = 2./(A2+A2)^ +7 2 ot . (S46) 



In order to apply these results for the Hamiltonians Eq. ( |S31[ ), (S39) describing the spectroscopy experiment, we 



make another transformation u<i = cxp [—i (uj cxt + kuj g ) ta z /2] into a frame rotating with a particular sideband k = 



—oo. ..oo. We can then use Eq. ( S44 ) for finding the steady-state population. One next makes the U2 transformation 
into the frame of each of the sidebands separately. If the peaks do not overlap, i.e., Slj <C w m , one can sum the 
occupancies due to each such operation: 



1 A 2 



where the Rabi frequencies arc given either by Eq. (S32) for the case of weak drive, or, more generally by Eq. ( S40 1 



In the latter case, the sideband Rabi frequencies are suppressed from those in Eq. (S32) due to the nested Bessel 
function dependence. 

Using these results, we can plot a prediction for the oscillatory Stark shift whose measurement data was shown 
Fig. 2d in the main text. We obtain an excellent agreement, as displayed in Fig. |S6[ The mirrored image is missing 
in the simulation since we have not included the offset quasiparticle in the simulation. The sideband visibility, 
proportional to its Rabi frequency, undulates with increasing n x according to the nested Bessel function dependence. 

Floquet calculation of the driven transmon 

According to the Floquet theorem^, the time-dependent Schrodinger equation, 

H{t)\^{t))=ihd t \*{t)), (S48) 
with a time-periodic Hamiltonian, H(t) = H(t +T), has a solution 

|*a(*)) = \u a {t))e- is ^ h , (S49) 
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FIG. S6: Analytical model of Fig. 2d in main text. The dashed line is the charge-qubit limit prediction, Eq. (S23l, 
however, with unrealistically low -Ej($dc/'l > o = 0.370) = 6.0 GHz. This is due to that the charge qubit model exaggerates the 
charge dispersion in the current transmon limit E.j > Ec- 



where \u a {t)) is periodic in time with period T — 2tt/uj. Substituting Eq. (S49| into Eq. (S48), we obtain the equation 
for the quasienergy state \u a (t)}, 



[H(t)-e a ]\u a (t)) =ihd t \u a {t)), 

which is just the time-dependent Schrodinger equation with the energy shifted by the quasienergy e c 
We expand the Hamiltonian and the quasienergy state with their Fourier components as 

#(<)=EE ft v ina; 'wi> 
w*))=EE c «e ima,t i«). 



(S50) 



(S51) 
(S52) 



Plugging these into Eq. ( S50 1 yields 



E E + <Wa,/jrofiwJ cf = ec n a . 

m /3 



(S53) 



Defining the Floquet matrix Hp as 



(an\H F \/3m) = h n a 







mhujS. 



n,m^a./3 : 



this eigenvalue problem can be written in a simple matrix form 

HpC = EC. 



(S54) 



(S55) 



Thus, we have transformed the finite-dimensional time-dependent problem to an infinite-dimensional time- independent 
problem. However, the number of Floquet blocks can be truncated depending on the strength of the periodic drive 
leading to computable eigenvalue problem. 

The Hamiltonian we consider is given in the Cooper pair basis by the qubit Hamiltonian (Eq. (S15l) with the gate 
charge driven sinusoidally, viz. Eq. (S14). This procedure allows for an exact treatment of the transmon qubit. The 
Fourier components of the Hamiltonian are 



= Hq , 

=4E c (n a -n)n x . 



(S56) 
(S57) 
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FIG. S7: Floquet results of the dynamical sweet spot, a, Transition frequency between the ground state and first excited 
state, as a function of gate charge modulation amplitude, b, Transition between the first and second excited states. c,d, 
Zoom-in of the regions marked in a,b, respectively. The family of curves depict changing no in steps of 0.05 between and 1/2 
as labelled. Notice the different scales in all panels. 



The Floquet Hamiltonian is thus given by 



H F = 



H [-i] 



(S58) 



where the dimension of and H^ 1 ^ is given by the number of Cooper pair number states used. In the modeling 
presented, 41 Cooper pair number states were used with 201 Floquet states. The size of the Floquet matrix was 
therefore 8241 x 8241. 



The transition frequencies plotted in Fig. 2d in main text and in Fig. S3 are then obtained as a difference of a 
proper pair of quasienergies^. We can also inspect the dynamical sweet spot (see also Fig. S5 1 in more detail, as 
displayed in Fig. [S7] For the lowest transition, the set of curves corresponding to different offset n cross at the same 
point at 0.2 % accuracy, whereas the next higher transition, the remaining dispersion is about 1 %. 

In order to obtain a numerically exact solution for the qubit population, one could in principle include the qubit 
excitation tone in the Floquet analysis as well, however, the straightforward diagonalization becomes intractable. 
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Driven response of the mechanical resonator 



With both dc and ac voltage applied to the mechanical resonator 

V g (t) = Vd c + Kc cos(oV), (S59) 



the driving force acting on it is 



where 



d_ fl 
dx \ 2 



C g (x)V^J =\F\cos(w g t), (S60) 



The induced motion is 



\F\ = V dc V ac ^- 
dx 



x(t) = , l f l cos K f + e ) = N M(Wjt + e) , (S6i) 



where 



tanO = "f' (S62) 

The gate charge due to the applied V g is 

n,(t) =^/- - "o + -^ — \x\ cosKt + 9) + -JL- cos( W ,t) (g63) 
=n + n x cos(uj g t + 0) + n v cos(uj g t) = n + r cos(uj g t + 9) . 
which has been expressed as a single sinusoid on the last line with 

r =\/nl + n x + 2n v n x cos 6 , (S64) 

and 

tantf = ^ sinQ (S65) 

n v + n x cos 

The phase shift of motion with respect to the driving voltage appears as an asymmetry in the Stark shifts in 
Fig. 2 in the main text. 

One also obtains a relation between the motional gate charge and phonon number N m — muj m x^ ms /h, 

N m = (^) n 2 x . (S66) 



The Stark shift, Eq. (2) in the main text, is at low phonon number 



Aw e _ f /27r^^7r 2 cos(27rno)(^) N m . (S67) 
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